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TERRESTRIAL ZONE DEBRIS DISK CANDIDATES IN H AND X PERSEI 
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ABSTRACT 

We analyze 8 sources with strong mid-infrared excesses in the 13 Myr-old double cluster h and x 
Persei. New optical spectra and broadband SEDs (0.36-8 /im) are consistent with cluster membership. 
We show that material with T^300-400 K and Ld/L* ^10 -10 -3 produces the excesses in these 
sources. Optically-thick blackbody disk models - including those with large inner holes - do not 
match the observed SEDs. The SEDs of optically-thin debris disks produced from terrestrial planet 
formation calculations match the observations well. Thus, some h and x Persei stars may have debris 
from terrestrial zone planet formation. 

Subject headings: planetary systems: formation planetary systems: protoplanetary disks 



1. INTRODUCTION 

Radiometric dating (Yin et al. 2002) suggests that 
the Earth had ~ 90% of its mass ~ 30-40 Myr after 
the formation of the Sun. Terrestrial planet formation 
models suggest that most of this mass accumulation was 
complete by ~ 10 Myr and was dominated by mergers of 
moon-mass oligarchs (Kenyon & Bromley 2006, hereafter 
KB06; Chambers 2001). How terrestrial planets form 
around other stars is unresolved. 

Optically-thin infrared (IR) emission from 
collisionally-produced debris is an important obser- 
vational signature of planet formation (Kenyon & 
Bromley 2004, hereafter KB04; Rieke et al. 2005, 
hereafter R05; KB06). For a 1-2 M Q star, terrestrial 
planet formation should produce debris (and hence 
strong 5-10 /im emission) for w 10-20 Myr. The 5-10 
/zm signature of debris production should then begin 
to fade (KB04; R05; KB06). While many debris disk 
surveys (R05; Gorlova et al. 2006; Su et al. 2006) 
have discovered evidence for Kuiper-belt analogues from 
24-70 fim excess emission around early A stars, there 
have been fewer searches for terrestrial zone debris 
disks around ~ 1-2 M© stars. Previous searches (e.g. 
Meyer et al. 2006; Silverstone et al. 2006) have also 
concentrated on protoplanetary disks < 5 Myr old, or on 
the evolution of classical debris disks (ages > 30 Myr) 
when planet formation is nearly complete and debris 
production is declining. 

The double cluster h and x Persei (d=2.34 kpc, ~ 13 
Myr old) provides an ideal probe of terrestrial planet for- 
mation. The clusters contain > 5000 stars with masses 
> 1.3 Mq, an order of magnitude larger population than 
other nearby, evolved clusters. The recent IRAC 3.6-8 
fim survey of h and x Persei by Currie et al. (2007; 
hereafter C07) discovered a large population of sources 
with IR excess emission in the IRAC bands (up to 4-8% 
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at 8j«m). C07 showed that the frequency of disk emission 
is larger for less massive stars and at greater distances 
from the central star. The IR excesses are bluer, K s - 
[IRAC]~ 0.5-1.5, than those for Class II T Tauri stars 
(e.g. Kenyon & Hartmann 1987; hereafter KH87), which 
suggests that the disks are less optically thick and per- 
haps in an evolved state. 

The goal of this paper is to constrain the source of 
the excess IR emission in h and x P er stars by modeling 
the SEDs and IRAC/MIPS colors of sources with high 
quality IR/optical photometry and spectra. In §2 we de- 
scribe our sample and explain our SED fitting. In §3 we 
show that the IR excess is best explained by debris from 
terrestrial planet formation. Both photospheric emission 
and optically-thick disk models fail to explain the SED 
slopes and IRAC colors. Optically-thick disk + inner 
hole models are marginally consistent with photometry 
< 5.8 \xm but are not consistent with either the 8/im 
([8]) data or MIPS data. 

2. INFRARED PHOTOMETRY AND SPECTRAL TYPING 

The C07 survey covered ~ 0.75 sq. deg. centered on h 
and x Persei. We select stars with J<16 and with 'strong 
excesses' at multiple IRAC bands, iT s -[5.8] > 0.5 (or K s - 
[4.5] > 0.5 in absence of [5.8] data) and K s -[8]>0.75. All 
of these sources have K s >13.5. In this sample, 1,915 have 
5cr detections at [5.8] (or [4.5]) and [8], lie along the h 
and x P er isochrone, and were not flagged for contamina- 
tion by background galaxies (see §2.3 of C07 for criteria). 
Eight of these sources (see Table 1) have 'strong excesses' 
and have high-quality optical photometry from Keller et 
al. (2001; UBVI) or Slesnick et al. (2002; V band only). 

We derive a worst-case probability that the [5.8] and 
[8] emission from these 8 sources are contaminated by 
z=0. 1-0.4 PAH-emission galaxies that escaped flagging 
and could masquerade as a 'strong excess' source. We 
assume that all sources from our population are as faint 
as our faintest source (K s =15) and that any superim- 
posed galaxy with [5.8] < 15 will trigger selection of one 
of the 1,915 sources. We also assume that all galax- 
ies are PAH-emission galaxies, that none were removed 
by flagging, and that each star covers an 'area' of 7r2 2 
square arcseconds. Using the galaxy number counts for 
the Bootes field from Fazio et al. (2004; ~ 1,165/sq. 
deg. for [5.8]<15), the expected number of contami- 
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nated sources is N p ~ 2.2 (tt2 2 x1165/(3600 2 )x1915). If 
we adopt the number density of PAH-emission galaxies 
for this field (Stern et al. 2005) instead of the number 
density of all galaxies, the predicted number of contam- 
inated sources is N p ^0.3. Thus, even under the most 
pessimistic assumptions no more than ~ 2/8 sources are 
contaminated. 

We obtained Hectospec (Fabricant et al. 2005) and 
FAST (Fabricant et al. 1998) spectra of these 8 sources 
on the 6.5m MMT and 1.5m Tillinghast telescope at F. L. 
Whipple Observatory during September-November 2006. 
For the FAST sources, we took ~ 10 minute exposures us- 
ing a 300 g mm -1 grating blazedat 4750 A and a 3" slit. 
These spectra cover 3800-7500 A at 6 A resolution. For 
each Hectospec source (the six faintest), we took three, 
10-minute exposures using the 270 g mm - * grating.^ This 
configuration yields spectra at 4000-9000A with 3A res- 
olution. The data were processed using standard FAST 
and Hectospec reduction pipelines (e.g. Fabricant et al. 
2005). The resulting spectra typically have moderate 
signal-to-noise, S/N > 30 per pixel (Figure la). 

To measure spectral types, we derived spectral indices 
of H a , H , tf 7 , H s , the G band (4305 i), and Mg I 

(5175 A) as in O'Connell (1973). We derived piecewise 
linear relationships between spectral indices and spectral 
type from the Jacoby et al. (1984) standards, similar to 
the method of Hernandez et al. (2004) . All sources have 
spectral types consistent with those expected for 13 Myr 
h & x Persei sources reddened by E(B-V)~ 0.52 (Bragg 
& Kenyon 2005) at a distance of 2.4 kpc. None are fore- 
ground M stars or background O/B stars. No sources 
have a H a equivalent width >10A indicative of gas ac- 
cretion. As an additional check on our spectral types, we 
fit the SED from U through K s to photospheric models 
in A vs. XF\ space. Spectral types derived from spec- 
troscopy and SED fitting agree to within 2-3 subclasses. 
To model their intrinsic SEDs through [8], the sources 
were dereddened in all bands. We use Bessell & Brett 
(1988) and Keller et al. (2001), respectively, to calibrate 
the V-I and U-B reddening from E(B-V), and the IR ex- 
tinction laws from Indebetouw (2005) to deredden the 
2MASS and IRAC bands. 

Finally, we acquired MIPS 24 (im photometry with an 
integration time of 80 seconds/pixel. The frames were 
processed using the MIPS Data Analysis Tool; PSF fit- 
ting in the IRAF/DAOPHOT package was used to ob- 
tain photometry adopting 7.3 Jy for the zero-point mag- 
nitude (Gordon et al. 2005). We derive the MIPS 5cr 
upper limits from the histogram of the number counts as 
a function of magnitude. The MIPS data are complete 
to ~ 10.5 and yield a 5<r detection for only one source, 
source 5, with [24]~9.9±0.05. We estimate the [24] red- 
dening from Mathis (1990): A [24] «0.025. Therefore, our 
one 5cr detection at [24] has if s -[24]~4.4. 

3. MODELING OF DISK EMISSION 

Figure 2 shows the source SEDs (filled diamonds). To 
compare them with stellar photospheres, we used op- 
tical/IR colors for each spectral type from Kenyon & 
Hartmann (1995; KH95) and converted from colors to 
XF\ in the Johnson-Cousins system. The zero-point flux 
from the models was scaled to match the dereddened V 
through J fluxes for each source (Allen 2001; Cohen et 



al. 2003). All sources have photospheric fluxes through 
[3.6]; one source (4) clearly has IR excess at [4.5]. All 
sources with [5.8] and [8] measurements show clear evi- 
dence for circumstellar emission. The SEDs are similar 
to some 'anemic' disks around A-F stars in IC 348 (Lada 
et al. 2006). 

To characterize the emission, we calculate the dust 
blackbody temperatures that match the observed [8] flux 
but do not overproduce the flux at other wavelengths. 
The range of dust blackbody temperatures (typically 
300-400 K) for each source imply terrestrial zone emis- 
sion. Derived disk luminosities are L^/D* ~ 10~ 4 -10~ 3 , 
comparable to debris disks (R05). 

To investigate the disk emission in more detail, we con- 
structed three disk models spanning a range of evolution- 
ary states for disks surrounding ~ 10 Myr old stars. The 
optically thick, non accreting, flat disk model (KH87) 
assumes that the disk passively reradiates photons from 
the star and extends from the magnctospheric truncation 
radius to > 2000 R*. In a second disk model, an 'inner 
hole' model, all disk material < 50 R± (~ 1 AU for R=2 
Rq, this was the best fit) from the star is removed from 
the optically thick, flat disk model. This model is similar 
to the morphology expected for 'transition' T Tauri disks 
(e.g KH95). 

Finally, we model the IR emission from an optically- 
thin debris disk. In a debris disk, modest IR excesses re- 
sult from small grains produced by destructive collisions 
between ~ 1 — 10 km planetesimals during the oligarchic 
growth stage of planet formation (see KB04; Kenyon & 
Bromley 2005). As in the KB04 calculations, the start- 
ing radial surface density (£) profile follows a minimum 
mass solar nebula, £ oc r~ 3 / 2 (Hayashi 1981). Plan- 
etesimal collisions are modeled in the terrestrial zone, ~ 
1.5-7.5 AU away from a 2M Q primary star. Disk emis- 
sion is tracked for ~ 100 Myr from the start of collisions. 
The resulting SED is computed and 2MASS and IRAC 
colors are determined. 

Figure 2 compares predictions from the three disk mod- 
els to the data (Figure 2). Optically-thin, debris disk 
models (dash-three dot line) consistently match the [5.8] 
and [81 IRAC fluxes. While the debris disk models can 



overpredict the flux at 
dictions at [5.8] and [8 



3.6] and [4.5] by - 25%, the pre- 
typically fall within 20% of ob- 
served values. In some cases (sources 2-4, 8) the agree- 
ment is excellent. Standard optically-thick disk models 
(long dashes, source 1 only) and the best-fit optically- 
thick disk + inner hole models (dashes) fail to match the 
SEDs. The optically-thick disk models overestimate the 
flux by > a factor of 10 beyond [3.6]. The inner hole 
model matches the observed SEDs through [3.6] but typ- 
ically overpredicts the [5.8] ([8]) fluxes by factors of 1.5 
(3). Furthermore, the optically-thick flat disk and inner 
hole models consistently predict a 24/ito flux ~ 10-20 x 
greater (or ~ 3 magnitudes redder) than the observed 
flux (one source) or the 5cr upper limits (5 sources). The 
debris disk model, which predicts 7f s -[24]~3, is consis- 
tent with the 24 fxm upper limits of the sources not de- 
tected with MIPS. The debris disk prediction is ~ 2.5 
times fainter (~ 1 magnitude bluer) than the observed 
if s -[24] color from source 5. However, planet formation 
calculations at 30-150 AU predict significant 24 /j,m emis- 
sion {K s -[24] « 3.5-4) by ~ 10-30 Myr around 2 M Q stars 
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(Currie et al. in prep., Kenyon & Bromley, in prep.), so 
regions beyond the terrestrial zone may contribute some 
[24] emission. More massive disks or those with flatter 
density profiles (e.g. r _1 ) than assumed here may also 
yield slightly larger [24] excesses. 

Figure 3 shows the results in K S -[A.5\/ K s -[8] and 
K s -[5.8]/K s -[8) color-color diagrams. In both figures, 
the IRAC colors of our sources (diamonds/triangles) lie 
halfway between the photospheric and inner hole mod- 
els (left square and dotted line). At the peak of debris 
disk emission (solid line) , the debris disk model also pro- 
duces IR colors about halfway between the photospheric 
and inner hole models. The colors show good agreement 
with debris disk predictions in i-T s -[5.8] and if s -[8], typi- 
cally within 0.2 mags, and fair agreement with the [4.5] 
data. Inner hole models for earlier (later) spectral types 
predict .Kg- [8] ~ 1.75 (3), or ~ 1-2 magnitudes too red. 
Many of the sources are, at best, marginally consistent 
with any star+background galaxy colors (dot-dashed en- 
closed region): 3/7 (4/7) are consistent with contami- 
nation in the K s -[4.S\/K s -[8] (K s -[5.8]/K s -[8\) diagrams. 
Only 2/8 are consistent with contamination in both plots. 
Even considering systematic errors/uncertainties (model 
uncertainties, photometric errors), the sources occupy a 
very small region of possible star+galaxy colors and are 
still only marginally consistent with contamination. 

The optically thick disk+inner hole model fails to 
model the IRAC/MIPS data in all cases; the debris disk 
model is consistent with the IRAC data and more consis- 
tent with the MIPS detection/upper limits. Thus, debris 
from planet formation is the most plausible explanation 
for near-to-mid IR emission from these sources. 

4. DISCUSSION 
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Near-to-mid IR-excess emission for our sample origi- 
nates from dust with T~300-400 K and Ld/L^ICT 4 - 
10~ 3 . Circumstellar debris produced from terrestrial 
planet formation is the most likely explanation for this 
emission. Future observations with the Spitzer Space 
Telescope and the James Webb Space Telescope can pro- 
vide stronger constraints on the SEDs of these terres- 
trial zone debris disk candidates and greatly expand the 
sample of ~ 13 Myr-old candidates in h and x Persei. 
Current IRAC/MIPS data require that the disks have 
little dust out to distances probed by [4.5] and are op- 
tically thin out to [24]. For a 13 Myr, KO (Fl) star 
with a luminosity of ~ 1.49 (7.4) L©, the disks are dust 
free out to ~ 0.23 (0.51) AU and optically thin to ~ 
6.6 (14.6) AU. Deeper MIPS [24] observations of these 
sources may more conclusively constrain the disk SED 
to compare with debris disk and transition disk models. 
Mid-IR spectroscopy may also aid source identification 
from the ~ 10-20/im flux levels and and the strength 
of the 10 fim silicate feature (see Sicilia-Aguilar et al. 
2007). 
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TABLE 1 

Sources analyzed in this paper (Dereddened, A^~1.62) 



# 


a 


8 


ST 


V 


U-B 


B-V 


V-I 


J 


H 


K s 


[3.6] 


[4.5] 


[5.8] 


[8] 


[24] 


1 


2:17:34.7 


56:47:25 


Fl 


14.15 








13.58 


13.51 


13.52 


13.38±0.01 


13.29±0.05 


12.97±0.05 


12.15±0.04 




2 


2:17:38.0 


56:53:51 


A9 


14.03 


0.09 


0.17 


0.06 


13.61 


13.55 


13.45 






12.55±0.06 


12.35±0.14 




3 


2:19:11.3 


57:09:53 


Fl 


14.72 


0.46 


0.37 


0.34 


13.88 


13.64 


13.64 


13.31±0.03 


13.23±0.02 


13.23±0.11 


12.89±0.08 


>10.5 


4 


2:22:06.4 


57:07:06 


F7 


15.44 






0.51 


14.62 


14.22 


14.28 




13.77±0.19 




13.38±0.12 


>10.5 


5 


2:19:07.7 


57:14:05 


F9 


15.77 


0.35 


0.63 


0.53 


14.65 


14.26 


14.29 


14.27±0.06 


14.19±0.09 


13.83±0.10 


13.03±0.09 


9.9 


6 


2:17:59.3 


56:53:45 


G9 


16.40 






0.45 


15.15 


14.60 


14.87 


14.71±0.07 


14.73±0.11 


14.27±0.13 


14.06±0.19 


>10.5 


7 


2:17:39.5 


57:20:26 


K0 


16.43 






0.76 


14.85 


14.43 


14.69 


14.43±0.06 


14.50±0.05 


14.29±0.18 


13.68±0.14 


>10.5 


8 


2:17:46.0 


57:10:18 


Kl 


16.72 






0.78 


15.12 


14.54 


14.43 


14.59±0.09 


14.35±0.10 


14.03±0.11 


13.51±0.12 


>10.5 
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Fig. 1. — FAST and Hectospec spectra for our sample, with spectral types ranging from Ff/A9 (top) to KO/Kf (bottom). The bottom 
source (8) has a 2000 count offset; spectra near the edges approach zero counts. 
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Fig. 2. — Spectral energy distributions of the sources. Filled diamonds denote detections; down arrows identify MIPS 5<r upper limits. 
The SED fits from V through [3.6] match photosphcric predictions (dot-dashed line). Sources have clear, strong IR excess at [5.8] and 
[8]. The dust blackbody temperatures indicate terrestrial zone emission. Both standard optically-thick fiat disk models (shown in the first 
plot only, large dashed line/open triangles) and optically-thick models with inner holes (dashed line/asterisks) overpredict the flux at [8]. 
Debris disk models (three dot-dashed line/open diamonds) appear to be consistent with the data. 




Fig. 3. — ii"s-[IRAC] colors for sources -triangles for sources earlier than G2, diamonds for sources later- and predictions for photosphcric, 
and standard optically-thick disk models (left and top-right squares), and a range of optically-thick disk + inner hole models (K0 to F0). 
The median errors in IRAC are shown in the lower left-hand corner. The photosphcric and optically-thick disk colors are for a spectral type 
of G2 and vary little with spectral type. Ovcrplottcd as a thick solid line are the expected colors for emission produced by optically-thin 
debris, a byproduct of terrestrial planet formation. The debris disk colors change due to the evolving dust production rate from planet 
formation. The dot-dashed region encloses possible colors for photospheric sources+background galaxy contamination; PAH galaxies 
have [5.8]-[8]=l-3 and any star+galaxy composite 'source' with [4.5,5.8]-[8]>l would have been flagged as a contaminated source and 
removed from our sample. The debris disk model fits to within ~ 0.3 mags in all cases; only 2/8 sources are consistent with star+PAH 
galaxy-contaminated colors. 



